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We develop a theory of collective spin dynamics triggered by ultrafast optical excitation of fer- 
romagnetic semiconductors. Using the density matrix equations of motion in the mean field ap- 
proximation and including magnetic anisotropy and hole spin dephasing effects, we predict the 
development of a light-induced magnetization tilt during ultra-short time intervals comparable to 
the pulse duration. This femtosecond dynamics in the coherent temporal regime is governed by 
the interband nonlinear optical polarizations and is followed by a second temporal regime governed 
by the magnetic anisotropy of the Fermi sea. We interpret our numerical results by deriving a 
Landau-Gilbert-like equation for the collective spin, which demonstrates an ultrafast correction 
to the magnetic anisotropy effective field due to second order coherent nonlinear optical processes. 
Using the Lindblad semigroup method, we also derive a contribution to the interband polarization 
dephasing determined by the Mn spin and the hole spin dephasing. Our predicted magnetiza- 
tion tilt and subsequent nonlinear dynamics due to the magnetic anisotropy can be controlled by 
varying the optical pulse intensity, duration, and helicity and can be observed with pump-probe 
magneto-optical spectroscopy. 
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I. INTRODUCTION 

The interaction between itinerant carrier spins and 
localized magnetic moments leads to carrier-mediated 
ferromagnetic order in a wide variety of systems, 
ranging from ferromagnetic semiconductors such as 
EuO, EuS, chrome spinels, or pyrochlorei to man- 
ganese oxides (manganites)^ and III-Mn-V ferromag- 
netic semiconductors^ili^. Such materials offer potential 
for novel spintronics applications^. Their magnetic and 
transport properties are intimately related and can be 
controlled by varying e.g. the carrier density, spin, and 
distribution. 

The non-equilibrium and dynamical properties of 
ferromagnetic semiconductors and magnetically-ordered 
systems are currently under investigation. For THz spin- 
tronic and magnetic devices, ultrafast information stor- 
age, recovery, and processing is required, e.g. the de- 
velopment of devices with sub-picosecond readout times 
of the magnetic states. This goal requires femtosec- 
ond spin manipulation and control. The physical pro- 
cesses that govern the magnetization dynamics during 
timescales shorter than the characteristic response times 
of the magnetic system are still under debate. Dur- 
ing such timescales, the validity of conventional ther- 
modynamic concepts for describing a magnetic system 
become questionable. Ultrafast pump-probe magneto- 
optical spectroscopy can shed light into this fundamen- 
tal problem. In these experiments, an ultra-short pump 
optical pulse excites optical polarizations, non-thermal 
populations, and carrier spins, which then trigger a mag- 
netization dynamics measured as function of time (Fara- 
day or Kerr rotation)^. 

To interpret such experiments, it is useful to distin- 
guish between different stages of time evolution of the 



photoexcited system. During the initial sub-picosecond 
regime, shorter than the dephasing times or optical pulse 
duration, the dynamics of the collective magnetization is 
triggered by optical polarizations and photoexcited car- 
rier spins. The response of the magnetic system is con- 
trolled by the nonlinear optical excitation and conven- 
tional thermodynamic concepts do not apply. In this 
initial regime, magnetic anisotropy can play a role by af- 
fecting the photoexcited carrier spin. After the pulse is 
gone, the light-induced quantum mechanical coherences 
decay and the photoexcited carriers relax by interact- 
ing with the hole Fermi sea. The carrier temperature is 
thus elevated above the lattice temperature within hun- 
dreds of femtoseconds. The hole population is eventu- 
ally described by a hot Fermi-Dirac distribution, which 
transfers its excess energy to the lattice within a few pi- 
coseconds. At the same time, the magnetic axes of the 
system change, due to e.g. the transient temperature 
elevation. Such quasi-thermal changes in the magnetic 
anisotropy are due to the Fermi sea carriers and should 
be contrasted to the magnetic anisotropy contribution 
of the non-thermal photoexcited carriers in the initial 
temporal regime. After the equilibration of the carrier 
and lattice systems, their common temperature relaxes 
via a slow (nanosecond) thermal diffusion process, which 
returns the magnetic system to its equilibrium configura- 
tion. In this paper we discuss a mechanism for coherent 
non-thermal magnetization manipulation and neglect all 
thermal effects due to elevated carrier and spin temper- 
atures. 

Most of the ultrafast magneto-optical experiments 
performed so far in magnetic metals, insulators, and 
semiconductors observed magnetization dynamics that 
could be interpreted in terms of light-induced time- 
dependent thermal effects. Following the observation of 
Ref§, many works focussed on ultrafast light-induced 
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demagnetization, which involves the time-dependent col- 
lective magnetization amplitude. The physical mecha- 
nisms that lead to quenching of the magnetization within 
a picosecond or less in materials ranging from transi- 
tion metals to III(Mn)V semiconductors are still un- 
der debate, but are mostly believed to be triggered 
by transient changes in the carrier and spin effective 
temperaturesi^'^iiSiiiiii^ii^'i^ However, transient mag- 
netic effects have also been observed in the initial non- 
thermal temporal regime, where the concept of car- 
rier temperature is not meaningful i^ii^ Light-induced 
changes in the magnetization orientation have also been 
observed in both metals and semiconductors and mostly 
attributed to transient changes in the magnetic easy 
axes due to temperature elevation.— Such 
thermal quasi-equilibrium effects induce a magnetization 
precession with a period ^lOOps. 

In most of the above experiments, the observed mag- 
netization changes result from the transient temperature 
rise following optical absorption and carrier relaxation. 
However, there is a limit in the magnetization speed 
that can arise from incoherent processes largely based 
on the heating of the magnetic system. Far more de- 
sirable is magnetization control based on coherent and 
non-thermal physical processes. The femtosecond coher- 
ent temporal regime offers the most flexibility for fast 
magnetization control limited only by the optical pulse 
duration. Coherent manipulation of the magnetization 
was demontrated experimentally in magnetic dielectrics 
and insulators22i2^iSi, while the interplay of coherent ex- 
citation and spin-orbit interaction in magnetic insulators 
was addressed theoretically in Refs. [25lf26l . 

Compared to other magnetically-ordered materials, 
IH-Mn-V ferromagnetic semiconductors offer certain ad- 
vantages and new possibilities for ultafast magnetization 
control and dynamics. These advantages stem from the 
carrier-induced nature of the magnetic order and the 
clear distinction between localized (Mn) and itinerant 
(valence band hole) spins. Static measurements have 
shown that III-Mn-V heterostructures are highly sensi- 
tive to external stimuli such as electrical gate, currents, 
or light i^iSSii^'^ A light-induced out-of-plane magneti- 
zation rotation towards the direction of propagation of 
a circularly polarized optical field perpendicular to the 
ground state magnetization was reported in Refs. Isills^ 
for Ga(Mn)As epilayers. Later ultrafast experiments in- 
terpreted their findings in terms of photoexcited carrier 
spin^ and thermali^ii^ effects (for a review see Ref. 0). 

Recently, Wang et.al.^'^ reported an enhancement of 
the magnetization amplitude and the ferromagnetic order 
in GaMnAs induced by the photoexcited hole population. 
This enhancement occurs on a ^lOOps time scale, fol- 
lowing the initial subpicosecond demagnetizationii and 
thermalization. More recently, Wang et.al^ reported 
the first observation of two temporal regimes of mag- 
netization dynamics in III-Mn-V semiconductors. The 
first regime lasts for a few hundreds of fs and is gov- 
erned by a quasi-instantaneous tilt of the collective mag- 



netization in response to optical excitation at high ener- 
gies (~ S.leV). In this femtosecond regime, a photoin- 
duced four-state ferromagnetic hysterisis was measured, 
which implies femtosecond detection of magnetic mem- 
ory states. This initial magnetization dynamics is clearly 
distinguished from the subsequent thermal regime, which 
is governed by magnetization precession on the lOOps 
timescale around the magnetic easy axes. 

The observations of RefH^ point out the need for a 
microscopic theory of collective spin dynamics in the ini- 
tial coherent regime of III(Mn)V semiconductors, which 
treats the nonlinear response of the magnetization to the 
ultrafast optical excitation. The theoretical prediction 
of a light -induced magnetic interaction resulting in a 
Kondo resonance in the nonlinear optical response and 
the pump-prob e sp ectra of doped semiconductors was 
reported in RefjSa. This light-induced many-body ef- 
fect is generated by a second order Raman-like process 
and should be most pronounced in the case of below- 
resonance photoexcitation in the transparency regime, 
where heating effects are suppressed. Ref. 37 suggested 
the possibility of inducing ferromagnetic order by excit- 
ing undoped paramagnetic II-Mn-VI semiconductors well 
below the optical absorption threshold. A brief descrip- 
tion of a microscopic mechanism for coherent ultrafast 
magnetization dynamics in III-Mn-V ferromagnetic semi- 
conductors was presented in Ref. [H, while Ref. [s^ ad- 
dressed the subsequent incoherent regime and attributed 
the ultrafast demagnetization to the scattering of the Mn 
spins with spin-flip excitations of the hot hole Fermi sea. 

In this paper we develop in detail a theory that de- 
scribes the ultrafast nonlinear response of the collective 
spins in ferromagnetic semiconductors. We present cal- 
culations for the most basic Hamiltonian that applies 
to a wide range of ferromagnetic semiconductor sys- 
tems. The magnetic exchange interaction and the cou- 
pling to the optical field are treated within the mean 
field approximation, which has been shown to describe 
well the ground state and thermodynamic properties 
of III-Mn-V semiconductors!^ as well as the ultarafast 
nonlinear optical response of semiconductors for strong 
photoexcitation.''"'*^ In Section |TT] we set up the prob- 
lem at hand, while in Section IIIII we derive the mean 
field equations of motion in the coherent limit in the case 
of several coupled valence bands. The dephasing effects 
due to the mixing of the hole spin states are described in 
Section [TVl with the Lindblad semigroup method^ This 
method allows us to treat consistently both hole spin re- 
laxation/dephasing and interband polarization dephas- 
ing resulting from hole spin-flip interactions. In Section 
|V]we extract a simple one-band model from the full the- 
ory and use it in Section lVTl to calculate the light-induced 
Mn spin trajectories and discuss the role of hole spin and 
interband polarization dephasing, magnetic anisotropy, 
and photoexcitation intensity and duration. Our nu- 
merical results predict an ultrafast tilt of the Mn spin 
away from its initial (equilibrium) value, which develops 
on a timescale comparable to the optical pulse duration. 
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The magnitude of this tilt is controlled by the optical 
pulse intensity and duration. The direction of the tilt 
and overall shape of the Mn spin trajectory depend sen- 
sitively on the hole spin dephasing and on the interband 
polarizations. In Section [VIII we interpret the above nu- 
merical results by deriving from the full theory an ef- 
fective Landau-Gilbert-like^^ equation of motion for the 
Mn spin using the adiabatic following approximation. We 
show that the effective magnetic anisotropy fields that 
govern the precession and relaxation within the above 
Landau-Gilbert picture acquire a time-dependent cor- 
rection determined by the optical pulse amplitude, the 
photoexcited interband polarizations, and the mixing of 
the hole spin states. This correction results in magne- 
tization dynamics on a timescale comparable to the op- 
tical pulse duration and may be interpeted in terms of 
light-induced precession and relaxation. After the pulse 
is gone, the Mn spin trajectory is controlled by the mag- 
netic anisotropics due to the Fermi sea (thermal) carri- 
ers. If the initial magnetization response to the optical 
excitation results in a sufficiently small tilt, the Mn spin 
precesses around its initial configuration, with a period 
determined by the zero-momentum magnon energy, as 
described by linearizing the equations of motion. On the 
other hand, for sufficiently high photoexcitation inten- 
sity, the initial light-induced magnetization tilt brings 
the magnetic system sufficiently far from the magnetic 
energy minimum so that nonlinear magnetic effects be- 
come important. In this case, the Mn spin evolution 
differs from a simple precession. Our theory predicts two 
distinct temporal regimes of magnetization evolution, the 
first of which is governed by the optical excitation while 
the second is governed by the magnetic anisotropics due 
to the thermal carriers. We discuss the signatures of the 
above transient and nonlinear magnetic effects in ultra- 
fast magneto-optical pump-probe spectroscopy and end 
with our conclusions in Section IVlIII The details of our 
calculations are presented in the five Appendices. 



II. PROBLEM SETUP 

We start with the Hamiltonian^i^ H{t) = H + HL{t), 
where H = + Kh + i?oxch is the Hamiltonian in the 
absence of optical excitation and {t) describes the cou- 
pling of the optical fields. 

Ke^Yl ^kneljkn , i^h = ^ e-k«^Lkn^-kn (1) 
kn kn 

describe the electron and hole band energies. The con- 
duction band electrons are created by the operator ej^^, 
where k is the momentum and n labels different conduc- 
tion bands. is the n-band dispersion (we set h=l). 
The valence holes are created by the operator where 
k is the hole momentum and n the band index. Their 
dispersion ej^^ is determined by the bandstructure. The 
Mn impurities act as acceptors, which create a hole Fermi 



sea in the valence band, and also provide randomly dis- 
tributed S — 5/2 spins, S^-, that are localized at positions 
Rj . These local moments interact with the hole spin via 
the antiferromagnetic Kondo-like exchange interaction 

-^cxch ^ ^ Sj • S^^/C ^ ^ ^^_kn^— k'n' 5 (2) 

where s is the hole spin operator, V is the volume, and 
(3 is the exchange constant. Finally, the coupling of the 
optical pulse is described by 

^Lit) = -J2 ^""'(04„^U„' + h.c, (3) 

nn'k 

where dnn'{t) — /i„„'£(i) is the Rabi energy, £{t) — 
£ exp [— t^/Tp] is the optical (pump) pulse, with dura- 
tion Tp, and iJ,nn' is the dipole transition matrix element 
between the valence band n' and the conduction band n. 
The optical pulse propagates along the growth direction 
z, which is perpendicular to the ground state spins. 

The ground state, thermodynamic, and transport 
properties of III(Mn)V semiconductors in the metallic 
regime (hole densities ~10^°cm~^) are well described 
by treating the magnetic exchange interaction, Eq.Q, 
within the mean field virtual crystal approximation!^ 
This approximation neglects spatial correlations and as- 
sumes uniformly distributed classical Mn spins, justified 
in the limit S —^ oo, where S is the Mn spin amplitude. 
The holes then experience an effective magnetic field pro- 
portional to the Mn spin. The valence band splits into 
two spin-polarized bands separated by the magnetic ex- 
change energy A = PcS, where c is the density of the 
localized Mn spins. The typical values of A observed in 
the GaMnAs and InMnAs ferromagnetic semiconductors 
that exhibit the highest critical temperatures are compa- 
rable to the Fermi energy, Ep ~100meV, of the ground 
state hole Fermi sea. 

In the ground state, the Mn spin points along the easy 
axis direction in order to minimize the total energy of the 
Fermi sea carriers: 

i5^(S) = ^i?^,nk„ (4) 

ki 

where i?^',-(S) are the eigensvalues of the Hamiltonian 
Kh + Hexch for given Mn spin S and nki are the popula- 
tions of the corresponding eigenstates. In III-Mn-V semi- 
conductors, this mean field total energy depends strongly 
on the orientation of S.''"''''^ This anisotropy is believed 
to mainly arise from the valence bandstructure, in par- 
ticular the spin-orbit coupling of different valence bands, 
which was taken into account in Refs,44,45 within the k-p 
envelope function approximation. The easy axis depends 
sensitively on the hole distribution among different va- 
lence bands. The calculated Eh{S) is well described by 
an expansion in terms of The following expansion 
describes the static experimental measurements:^^ 

Efi = Kc{S'^Sy + S'^S'^ + SyS'^) + KuS'^ + KuzS'^, (5) 
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where S = S/S, Kc is the lowest cubic anisotropy 
constant, and and K^z are the first order uniax- 
ial anisotropy constants, attributed to strain, whose ori- 
gin is still being debated.''^ For — K^z — 0, the 
ground state magnetization is either parallel to the axes 
±a;, iby or ±2;, if Kc > 0, or points along the diagonals 
= \y\ = l-^li if < 0. For the parameters of inter- 
est in Ga(Mn)As epilayers, < 0, and the easy axis 
direction is determined by the magnitude and sign of the 
uniaxial anisotropy constants K^^ and K^^^z- A sufficiently 
large Ky^z > ensures that the ground state magneti- 
zation lies within the x-y plane as observed experimen- 
tally. The observed temperature and hole concentration 
dependence of the easy axis suggests that Ku > 0^^ For 
Ku > \Kc\, the easy axis points along the ±y axis, while 
for Ku < \Kc\ it points at an angle (j) from the x-axis, 
where cos 20 — Ku/\Kc\. 

Even though k • p and mean field theory explain the 
main anisotropy effects, the interpretation of the ob- 
served changes of the in-plane easy axis with temper- 
ature and hole concentration require further theoretical 
investigations. A complete microscopic theory that in- 
cludes the magnetic anisotropy becomes even more com- 
plicated in the case of ultrafast optical excitation. As 
observed experimentally,— the magnetization responds 
to the photoexcitation of high energy (~3.1eV) valence 
bands states, which are empty in the ground state and 
lie far from the Brillouin zone center, well before carrier 
or spin thermalization. In this initial highly nonequi- 
librium regime, a fully microscopic theory must address 
both the valence bandstructure at ~3eV energies and the 
coherent and non-thermal effects. Here we address the 
general magnetic semiconductor system and include the 
magnetic anisotropy effects by using the Fermi sea energy 
expression Eq.©. The time dependent response of the 
Mn spin S to the nonlinear optical excitation, discussed 
in the next section, changes Eh{S), which leads to an 
additional complex nonlinear magnetization dynamics. 



III. EQUATIONS OF MOTION 

To describe the ultrafast optical and spin response, we 
proceed in two steps. In this section we derive the mean 
field equations of motion, while in the next section we 
derive the dephasing contributions due to the hole spin- 
flip interactions. The components of the density matrix 
(p) are obtained from the equations of motion 



^^t{p)= {[p,H{t)])\ 



HF 



^dt {p)\ 



relax ' 



(6) 



where the last term describes the dephasing and relax- 
ation contributions. By factorizing all higher density ma- 
trices (Hartree-Fock approximation), which couple due 
to the many-body exchange interaction, we obtain a 
closed system of equations for the optical polarizations, 
spins, and carrier populations and coherences. 

We are interested in calculating the macroscopic mag- 
netization measured in ultrafast pump-probe magneto- 



optical experiments. This is dominated by the contribu- 
tion of the average Mn spin 



(7) 



From Eq. ^ we obtain within the mean field approxima- 
tion the equation of motion that describes the magneti- 
zation dynamics: 



as = S X 



H 



(8) 



where H is a magnetic field. The right hand side (rhs) 
of the above equation describes the precession of the Mn 
spin around an effective time-dependent magnetic field 
determined by H and the mean hole spin 



(9) 



where s„„' are the matrix elements of the hole spin op- 
erator s between the valence band eigenstates. Eq.® 
corresponds to the Landau-Gilbert picture of magneti- 
zation dynamics*'^ and conserves the amplitude of S. 

The hole spin, populations, and inter-valence band 
coherences are described by the density matrices 
{h}_y^^h-\^n') , whose equations of motion read 

m' 

kn^km' / 



(10) 



The second line on the rhs of the above equation de- 
scribes the change in the hole states due to the exchange 
interaction. The third line describes the excitation of 
hole coherences and populations by the optical pulse, via 
second-order Raman processes. Similarly, we obtain for 
the electron populations and coherences 

«9t(e[„ek„') = (Ckn' - (eLek«') 



■i 9t(e[„ek„' 



(11) 



The total carrier populations with given momentum k, 

K = E<4„ekn) , iVk' - E<^-kn/^-kn>, (12) 
n n 

are not affected by the magnetic exchange interaction, 
which within the mean field approximation only changes 
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the spin. They satisfy equations of motion similar to the 
Bloch equations: 

nn' 

(13) 

and 

dtN^ = d;,„(t)(/i_k„ek„'> - K/^f, (14) 

nn' 

where T/* and Tf are the spin-independent population 
relaxation times, due to the carrier-carrier and carrier- 
phonon scattering processes which thermalize the carrier 
system. The hole population relaxes to the thermal dis- 
tribution 

Using Eq.® for the total hole spin, we obtain from 
Eq. after some algebra that 

Stsij = PcS X sjj + /,7ihk(<) + i{[Kh,s'^]) + dts'^l^^^^ . 

(15) 

The first term on the rhs describes the precession of the 
hole spin around the mean field created by the Mn spin. 
The second term describes the photoexcitation of hole 
spin via second-order Raman processes, which are deter- 
mined by the transition matrix elements (selection rules) 
and the interband optical polarizations: 

hk(0 = 2 ^ d*„^,{t)Sm'n'{h-kn'eWn)- (16) 
nn'm' 

The third term on the rhs of Ea. (|15p arises when the 
valence band states are not eigenstates of the hole spin: 
[sjj, Kh] ^ 0. In this case, the hole spin dynamics de- 
pends on the valence band energies and is determined 
by the individual coherences, Eg. (fTO|) . between the dif- 
ferent valence bands. In III-Mn-V semiconductors, the 
hole spin is not conserved due to the strong spin-orbit 
interaction. Recalling Eq.®, we conclude that the mix- 
ing of the hole spin states and the valence bandstructure 
modifies the effective magnetic field experienced by the 
Mn spin as described by the third term on the rhs of 
Eq.dlll). 

The photoexcitation of carrier spin, which triggers the 
magnetization dynamics, is governed by the interband 
optical polarizations 

Pwnn' = {h-knekn') , (17) 

determined by the equations of motion 

idtPkrin' = {£kn' + ^kn ^ " 'i/T2) Pknn' 

-|-/5cS • ^ ^ ^nmPk.mn' 
m 

-dn'n{t) [l - (ALkn^-kn) " (eL'^kn') 

m^n m^n' 
+idt Pknn'lelax ■ (18) 



In the above equation, ujp is the pump optical pulse cen- 
tral frequency and T2 describes the spin-independent 
contribution to the polarization dephasing. T2 can be 
quite short due to the disorder, which relaxes the momen- 
tum conservation by introducing a one-body potential. 
For weak disorder, this effect can be treated by introduc- 
ing a dephasing time comparable to the momentum scat- 
tering time. The second line on the rhs of the above equa- 
tion is due to the change of the hole states induced by the 
effective magnetic field f3cS{t). The third line describes 
the Pauli-blocking nonlinearity (Phase Space Filling)"*-, 
while the fourth line describes the contribution of carrier 
coherences between the different bands. Finally, the last 
line describes the spin-dependent polarization dephasing 
due to the hole spin-flip interactions, discussed in the 
next section. We note from the above equations that the 
interband polarizations, valence band coherences, and 
hole states depend on the effective magnetic field propor- 
tional to S{t). This spin can deviate significantly from 
its ground state configuration in the case of strong pho- 
toexcitation, which in turn changes the hole states as 
compared to the ground state. Such light-induced de- 
viations from equilibrium increase with photoexcitation 
intensity and are described non-perturbatively by solv- 
ing numerically the above system of coupled equations 
without expanding in terms of the optical field. 



IV. HOLE SPIN AND POLARIZATION 
DEPHASING 

A. Lindblad formalism 

The equations derived in the previous section do not 
include dephasing and relaxation contributions. Within 
the semiconductor Bloch Equations)^ such effects are 
treated phenomenologically by introducing effective de- 
phasing and relaxation times to the polarization and pop- 
ulation equations of motion. Here we must also consider 
the dynamics of the carrier and Mn spins, which within 
the mean field approximation precess around each other 
while conserving their magnitudes. In the III-Mn-V sys- 
tem, the hole spin relaxation is strong, mainly due to the 
spin-orbit interaction. For example, the combination of 
spin-orbit coupling in the valence band and disorder- 
induced scattering between momentum states, leads to 
hole spin dephasing and relaxation times of the order of 
lO's of fs42ii^ Such times are comparable to the hole spin 
precession period around the Mn spin and cannot be ne- 
glected. 

Carrier spin relaxation is often described within the 
spin Bloch equations.^^'^*' However, hole spin relaxation 
also leads to interband polarization dephasing, and both 
effects must be treated on equal footing for our pur- 
poses here. For this we use the Lindblad semigroup 
description of dissipative quantum dynamics. Under 
the general assumptions of linear coupling between bath 
and system operators and Mar kovian/relaxation time ap- 
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proximation, as well as density matrix positivity and 
semigroup-type time evolution, the relaxation contribu- 
tion to the density matrix equation of motion can be 
expressed in the form^ 



dtp\re 



{Ll„,,Li,'m'P) - (p4'™'^k'™')l ,(19) 



where r± is the spin dephasing rate. The Lindblad op- 
erators £km must be chosen to describe the relaxation 
processes at hand. In our case, the hole spin must relax 
towards the direction antiparallel to the Mn-spin S{t). 
As we show below, this can be achieved by choosing 



(20) 



where h^^^ ( ^km^) creates a hole with spin parallel 
(anti-parallel) to the Mn spin S{t) and the index m labels 
the different basis states. 



B. Hole spin dephasing and relaxation 

First we derive the hole spin dephasing and relaxation 
described by the Lindblad operator Eg. ([^0)1 . Noting that 
the hole spin relaxes to a direction antiparallel to the Mn 
spin, it is useful to introduce its components parallel and 
perpendicular to the unit vector S: 



the basis of hole spin eigenstates discussed above and 
define the interband electron-hole amplitudes 



kmncT 



(26) 



where a =1,1- Using Egs. fTO)) and (^0]) . we derive 
in Appendix (|B]) the following expression for the spin- 
dependent polarization dephasing: 



dtPkmnI \relax 



-r, X 



Pi 



kmnf 



2 ^ 



— (1-^1 
2 ^ ' 



h N 

km/ 



*(^ X sjj„)j 



Pi 



kmnX 



(1 - NtJ + ^iS X stJ- 



(27) 



dtPlimnilrelax ~ ^ 



kmn]. 



2 ^ 



Pi 



kmni 



(28) 



where, for any vector A, A± = i lAy. The above 
dephasing contribution depends nonlinearly on the pho- 
toexcitation, via the photoexcited hole contribution to 
the population and spin s'' and the light-induced 
changes in the Mn spin S. These nonlinearities corre- 
spond to excitation-induced dephasing induced by hole 
spin-flip correlations. 



where 

krj 



*krr 



S ■ sjj„ and 



Ss 



kmll ' 



(21) 



s;:„^ = S X (sU„ X Sj . (22) 

As derived in Appendix [A} the Lindblad operator Eq. ((20|) 
gives the following spin-Bloch equations: 



^limA.\relax 



Fis 



km_L' 



'km I , 



where Fn = 2F^ is the spin relaxation rate and 



(23) 
(24) 

(25) 



is the quasi-equilibrium hole spin value, corresponds 
to the maximum spin, Smax, of TV'' holes, given by the 



relation Sn 



C. Spin— dependent Polarization Dephasing 

The hole spin dephasing also dephases the interband 
optical polarizations. To describe this, it is useful to use 



V. SINGLE BAND APPROXIMATION 

The general theory derived in the previous sections 
may be used to treat the bandstructure relevant to the 
particular ferromagnetic material of interest. Below we 
discuss the general features of the dynamics, their phys- 
ical origin, and their sensitivity to different parameters 
by extracting from the general theory a simplified one- 
band model that captures the essential physics common 
in all materials. We therefore neglect bandstructure par- 
ticularities, such as the nature of the high energy states 
far from the Brillouin zone center, excited by the ~3.1eV 
pump of Refi35„ or impurity bands. We assume that 
the magnetization dynamics is mainly triggered by opti- 
cal transitions between a single valence and conduction 
band of spin-f and spin-| states, whose mixing is de- 
scribed phenomenologically with the Linbdlad approach. 
In the basis of carrier spin eigenstates and assuming a 
single band, Eq.([3]) reduces to 



Hr 



E 



''ki''-kT 



^'-E^kT^-ki + 'i-C- (29) 



where d±{t) — p±£{t). The Rabi energy d+ describes 
the coupling of the right-circularly polarized component 
of the optical field, while c?_ describes the coupling of the 
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left-circular ly polarized component. The corresponding 
interband transition matrix elements /i± depend on the 
bandstructure and the admixture of spin-f and spin- J. in 
the band states that mostly contribute to Eq.Q for the 
energies of interest. Denoting 



P+ 



(/i-k<Tekj,),Pk£7 = (/^-kCTekt)- 



(30) 



and using Egs. f??]) and ([^5]) and the results of Section 
mil we obtain the following coupled equations of motion 
for the interband optical polarizations: 



idtP+ 



kT 



-■^ki 



= -d+{t) 



- Ak+P, 



kT 



idtP^^. - (r!k + Ak^)PkT - ^k-Pk^ 

^d.{t) 4_ + d+{t) s^_, 



where 



1 



^k|| 



(31) 
(32) 
(33) 
(34) 

(35) 



gives the Mn spin-independent and 



Ak = -[/3c5-zri(l-0] +riSxs(j (36) 

the Mn spin-dependent contribution to the e-h pair en- 
ergy and dephasing. Ak also determines the coupling 
between the interband polarizations, due to the mix- 
ing of the hole spins by the magnetic exchange inter- 
action and the dephasing. The above e-h pair energies 
depend on the light -induced deviations of the Mn and 
hole spins from their ground state configurations and 
on the photoexcited hole populations, which give non- 
linear contributions to the optical polarization. The rhs 
of Eas. (PT|) and ([55]) describes the Pauli-blocking nonlin- 
earities (Phase Space Filling)^, while the rhs of Eqs. ([5^ 
and (j34p describes the contribution of carrier spin coher- 
ences, ground state or photoexcited. The above polariza- 
tion equations of motion treat the effects of the mixing 
of the spin-| and spin-| states by spin-orbit or other 
spin-flip interactions by using Eqs. ipT)) and 

Finally we turn to the equations of motion for the 
carrier spins. The hole spin dynamics is described by 
Eq. (|15|) . Within the one-band approximation, the mix- 
ing of the spin-f and spin- j hole states leads to the spin 
dephasing derived in Appendix |^ The second term on 
the rhs of Eq. (fT5|) describes the photoexcitation of hole 
spin as determined by the interband polarizations and 
hk(i), Eq. fTO]) . In the basis of hole spin eigestates con- 



sidered in this section, 

hk.(i) = rf;(t)Pk+(<) + rf*_(i)Pk^(i) 



hky(0 = -i 



MP+^{t) - d*_it)p-it) 



hk4i) = dl{t)P+{t)-d*Jt)P^^{t). (37) 

The above equation describes a second-order nonlinear 
optical process where the excitation of a spin-| hole- 
spin-| electron pair is followed by the de-excitation of 
a spin-| hole-spin-| electron pair. This second-order 
(Raman) process is induced by the right-circularly po- 
larized component of the optical field, while the left- 
circularly polarized component induces an analogous pro- 
cess involving the electron spin-| states. The above pro- 
cess requires a nonzero polarization Pkj- This is possi- 
ble on the one hand due to the coupling of the spin-f 
and spin- J, hole states (and hence P^ with Pj^, third 
term on the Ihs of Eq. ([5^ ) described by Ak, due to the 
magnetic exchange interaction and the hole spin dephas- 
ing processes (e.g. spin-orbit), and on the other hand 
due to the presence of spin-f-spin-J, hole spin coher- 
ence in the ground state. Finally, Ea. (P7|) describes the 
photoexcitation of spin-polarized hole populations and 
a hole spin z-component via the excitation and subse- 
quent de-excitation of a spin-| hole and spin-| electron 
pair (right-circularly polarized component) or a spin-J, 
hole and spin-| electron pair (left-circularly polarized 
component). The above coherent nonlinear effects oc- 
cur during the optical pulse and do not require the ab- 
sorption of light. They can also be induced by photoex- 
citing the system below resonance, in the transparency 
regime, provided that the photoexcitation intensity is 
sufficiently high to achieve an observable effect. Such 
below-resonance photoexcitation is advantageous since 
undesirable effects such as transient heating can be sup- 
pressed and thus the speed of a possible device can be 
maximized. The equation of motion for the electron spin 
has a form similar to Eq. p^ : 



(dt + r^Jsl^ ImKit) (38) 

where is the electron spin dephasing rate and de- 
scribes the photoexcitation of conduction electron spin: 

hL(^) = d*_it)P+it)+d*,it)P^^{t), 
Kyit) = ^[dl{t)P^^it)-d*_it)P+it)^ 

KM - -Kit)P+it)^d*_it)p-^it)]. (39) 

In the next section we calculate the Mn spin dynamics 
by solving the above system of coupled equations non- 
perturbatively, which allows us to treat large deviations 
of the Mn spin from its ground state configuration in- 
duced by strong photoexcitation. 



VI. NUMERICAL RESULTS 

Within the mean field approximation, light-induced 
Mn spin dynamics is triggered initially by the photoex- 
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(a) 




FIG. 1: (Color online) Mn spin trajectories without mag- 
netic anisotropy in the three regimes of spin and polariza- 
tion dephasing discussed in the text for three pulse dura- 
tions Tp and two different Rabi energies: (a) d+{0)—20meV , 
(b) d+(0)=60meV. Thick circled line: Tp^lOOfs, Dashed line: 
Tp=250fs, Thin solid line: Tp=500fs. 



citation of a hole spin component perpendicular to the 
ground state Mn spin. In the simple one-band approxi- 
mation, the magnitude of such a spin component is deter- 
mined by the ratio d+ / d- between the Rabi energies that 
describe the coupling of the right- and left-circularly po- 
larized components of the optical field. This ratio is de- 
termined by the selection rules and the nature of the 
bands that contribute to the magnetization dynamics. 
Within our simplified one-band model, /c?_ also re- 
flects the magnetic anisotropy of the system due to the 
bandstructure. In the absence of a complete theoretical 
understanding of such bandstructure effects, especially 
for the high energy (^S.leV) transitions observed ex- 
perimentally to trigger ultrafast magnetization dynam- 
ics in Ga(Mn)Asr^ we consider here the extreme case 
of d- = 0, which corresponds to right-circularly polar- 
ized light. For all the calculations presented in this pa- 



per, we consider a magnetic exchange energy /3cS'=125 
meV comparable to the Fermi energy Ep=100 meV, frac- 
tion of initial holes 0.33 of the Mn impurities, a hole 
mass rrih = l.lbrrie, and a carrier thermalization time 
Ti=165fs. Our results are not very sensitive to the above 
parameters, with the exception of the exchange inter- 
action (3 that changes the magnitude of the effect. To 
make the connection between the calculated quantities 
and the experiment, we note that ultrafast magneto- 
optical pump-probe spectroscopy can be used to deduce 
the time evolution of the z-component of the magnetiza- 
tion, - S,{t)J.^ 

We start with our results in the absence of magnetic 
anisotropy. Fig [1] shows the light-induced time evolution 
of the Mn spin and its dependence on the hole spin and 
polarization dephasing, the Rabi energy, and the optical 
pulse duration. By comparing three different dephasing 
regimes in Fig[Tl one can see that the shape of the Mn 
spin trajectory depends sensitively on the hole spin and 
polarization dephasing. Regime I corresponds to very 
short dephasing times, l/r^=21fs and T2=10.5fs, which 
are typical in III(Mn)V semiconductors.^* Regime 2 cor- 
responds to r2=330fs with very short hole spin dephasing 
l/r_L=21fs. Finally, Regime 3 neglects the spin dephas- 
ing altogether, Fj, = 0, T2=330fs, which assumes that 
the hole spin is a good quantum number. As can be seen 
in FiglU the dephasing changes qualitatively the shape of 
the Mn spin trajectory, from the precession within the x- 
y plane expected in the absence of magnetic anisotropy 
to a complex magnetization tilt out of the x-y plane, 
in the direction of pulse propagation. The magnitude 
of this out-of-plane tilt and the Mn spin z-component 
are enhanced by increasing Fj^ or the pulse duration Tp 
and by the short T2. Furthermore, the Mn spin compo- 
nent within the x-y plane rotates in a clockwise direction, 
opposite to the counter-clockwise direction expected for 
an effective magnetic field along the +z axis of pulse 
propagation. For right-circular-polarization, the latter 
is the expected direction of the photoexcited hole spin 
in the absence of magnetic exchange interaction or spin 
dephasing. However, our numerical results show that the 
photoexcited hole spin z-component is in fact negative, 
due to the magnetic exchange interaction and the mix- 
ing of the hole spin states. The dependence of the tra- 
jectories in Fig[l] on the pulse duration shows that the 
light-induced Mn spin dynamics terminates soon after 
the end the photoexcitation, as expected in the absence 
of magnetic anisotropy. Furthermore, the magnitude of 
the Mn spin tilt increases with photoexcitation intensity. 
The spin dynamics of Fig [1] describes the fundamental re- 
sponse of the spin system to the optical excitation within 
the mean field approximation. 

The Mn spin response is seen more clearly in Figl^l 
which shows the development with time of the y- and z- 
components of S. Both of these components vanish in the 
ground state. They develop on a timescale determined 
by the pulse duration and have comparable magnitudes 
for sufhciently large Fj^. The importance of the hole spin 
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FIG. 2: (Color online) Mn spin components without mag- 
netic anisotropy for pulse durations and Rabi energies as in 
FiglU l/ri=21fs, r2 = 10.5fs. Dashed line: Sy{t). Dotted 
line: S^{t) 

dephasing can be seen by comparing Figl^l obtained for 
very short spin and polarization dephasing times, with 
FigEl obtained for Fj^ = 0. The most striking difference 
is the very small magnitude of the out-of-plane compo- 
nent Sz {t) when the hole spin is conserved (Fj^ = 0) . The 
hole spin follows the overall Mn spin and remains more or 
less antiparallel to S{t) at all times, with the exception of 
a small component perpendicular to S that triggers the 
spin dynamics. 

We now turn to the effects of the magnetic anisotropy 
on the spin dynamics. As discussed above, the magnetic 
anisotropy of the thermal carriers leads to preferred Mn 
spin directions that minimize the total Fermi sea energy 
Eh{S), Eq.([5l). Its effects on the Mn spin dynamics can 
be treated phenomenologically by adding the magnetic 
field 

as 

to Eq.®.''^"'^^ Such an anisotropy field causes a nonlinear 
rotation and relaxation of the Mn spin until it aligns with 
the easy axis direction that minimizes Eh{S), so that 
H = 0. We would like to distinguish between the above 
effect, due to the thermal carriers, and the anisotropy 
effects on the photoexcitation. The latter can be treated 
microscopically as discussed in Section III and determine 
the photoexcited hole spin (see e.g. Eq. irSl) ). 

In Ga(Mn)As, Kc < 0, X„ > 0, and \Kc\ > K^-^^'^'^ In 



FIG. 3: (Color online) Time evolution of Mn spin components 
for the parameters of Fig[2] except for rx=0, r2=330fs. 



this parameter regime, there are two degenerate ground 
states, which correspond to Mn spin pointing at angles 
(f) from the x-axis such that cos 20 — Ku/\Kc\- The 
system can be prepared so that, prior to the photoexci- 
tation, the Mn spin points along either one of the above 
two easy axes. In the absence of magnetic anisotropy, 
both these initial conditions would result in the same Mn 
spin trajectory. FigsHI^a) andlH^b) compare the three Mn 
spin components as function of time, starting with Mn 
spin pointing along the two different ground state con- 
figurations. To make this comparison more meaningful, 
in both cases we chose the x-axis of the coordinate sys- 
tem to coincide with the initial magnetization direction. 
FigHl clearly shows that, as the Rabi energy increases, 
the two initial conditions lead to different time evolu- 
tion, which implies that the two magnetic ground states 
can be distinguished. As can be seen by comparing FigU 
with FiglU the Mn spin dynamics depends on the hole 
spin dephasing, especially as the photoexcitation inten- 
sity increases. 

The dependence of the Mn spin dynamics on the ini- 
tial magnetic state becomes more clear in FigsIB] and [71 
which plot the Mn spin trajectories for the two initial 
conditions in the same coordinate system, whose x-axis 
coincides with the initial magnetizations. Two different 
temporal regimes can be clearly distinguished. The first 
regime lasts for a time interval comparable to the pulse 
duration. Here the dynamics is governed by the opti- 
cal pulse intensity, duration, and helicity, as well as by 
the ratio /i+ and the mixing of the hole spin states 
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FIG. 4: (Color online) Magnetic anisotropy effects on the Mn 
spin components for strong dephasing l/rx=21fs, r2=10.5fs, 
pulse duration tp—500is, and two Rabi energies. Solid line; 
Sx{t). Dashed line: Sy{t). Dotted line: Sz{t). (a) and 
(b) compare the time evolution for initial condition along 
the two easy axes, taken as the x-axis. -ftrc=-0.0175meV, 
7i'„^=0.0252meV, 74:„=0.010meV. 



that depend on the magnetic anisotropy. On the other 
hand, the magnetic anisotropy field Eq.([3D]), due to the 
thermal carriers, plays a very minor role in this initial 
temporal regime for sufhciently large Rabi energies that 
exceed the anisotropy constants. The trajectories for the 
different initial conditions coincide in this sub-picosecond 
regime, which gives an ultrafast magnetization tilt deter- 
mined by the response to the hole spin photoexcited by 
the coherent spin Raman processes described by Eqs. p^ 
and ((57|) . As can be seen by comparing Figs (Hand [71 the 
shape of the Mn spin trajectories is strongly influenced 
by the spin and polarization dephasing. 

FigsE] and [7] show a rather abrupt change in the Mn 
spin trajectory at times i ~ r^. At later times, the tem- 
poral evolution is determined by the magnetic anisotropy 
field Eq. pO)) . The shape of the trajectory during many 
picoseconds depends critically on the magnetization tilt 
that develops during the femtosecond initial stage. If 
this tilt is sufficiently small, the Mn spin dynamics is 
described by linearized equations of motion, expanded 
around the equilibrium spin values that minimize the 
Fermi sea magnetic energy Eh(S). This harmonic os- 
cillation corresponds to zero-momentum magnon excita- 
tions, whose frequency is the same for the two easy axes. 
Since the magnitude of the magnetization tilt can be con- 
trolled by the optical pulse intensity and duration, we 
conclude that for sufficiently short pulses and sufficiently 
small Rabi energies, rf+Tp < 1, the Mn spin trajectories 
are very similar for both initial conditions. They corre- 
spond to magnetization precession around the easy axis 
(magnons). However, with increasing photoexcitation in- 
tensity and duration, the deviation from the easy axis due 
to the sub-picosecond magnetization tilt increases, and 



FIG. 5: (Color online) Time evolution of Mn spin components 
for the parameters of Fig|4] except for rx=0, r2=330fs. 

eventually the Mn spin dynamics cannot be described by 
expanding i?/i(S) around its minimum. Such nonlinear- 
ities result in a complex trajectory, determined by the 
full nonlinear equations Eqs.lH]) and (|40)) . which can dif- 
fer substantially from a simple magnon precession. As 
demsonstrated by Figsin]and[71 in this nonlinear dynam- 
ics regime, the Mn spin trajectories depend on the initial 
condition. One can then distinguish between the differ- 
ent magnetic states, which can provide the basis for an 
ultrafast magnetic memory read-out scheme with speed 
limited only by the optical pulse duration. Such a scheme 
would be based on controlling the ultrafast response of 
the spin system via the optical pulse intensity, duration, 
and helicity. Furthermore, the strong dependence of the 
trajectory shape on the spin dephasing implies that the 
valence bandstructure and mixing of the hole spin states 
by the spin orbit interaction plays an important role. In 
the next section we interpret our numerical results by de- 
riving from the full theory an effective Landau-Gilbert- 
like equation for the Mn spin after expanding around the 
adiabatic limit. 



VII. INTERPRETATION: LIGHT-INDUCED 
RELAXATION AND RE-ORIENTATION 

The numerical solution of the full mean field equations, 
discussed in the previous section, shows that the hole spin 
follows the Mn spin more or less adiabatically. This is due 
to the much faster hole spin precession and relaxation as 
compared to the time-dependent changes in S{t). To 
interpret the full numerical results, we therefore expand 
around the adiabatic limit. The derivation presented in 
Appendix [U] then gives a Landau-Gilbert-like^ equation: 

dtS (1 - 7s|f)S X H + as|S x (^S x 

- aS X h{t) - 7S X (s X h(t)) , (41) 
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FIG. 6: (Color online) Mn spin trajectories with initial con- 
dition along the two different easy axes (I and II) for the 
parameters of Fig|4] (a): d+=40meV, (b): d+=60meV. 
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FIG. 7: (Color online) Mn spin trajectories for the same pa- 
rameters as in Fig|6]except for rx=0, T2=330fs. 



where the effective magnetic fields h (Ea. ([57)) ). 

hW^-^E ^"^hk(0, (42) 

k 

and H, Eq.([3D]), are due to the optical excitation and 
magnetic anisotropy respectively. In the above equation, 

P ((3cS + f3s'i + Hn) 

7 = ,2 ^ (43) 

i^(3cS + Ps\^ + H^\j +Tl 

where _ff|| = H ■ S, while 



The first two terms on the rhs of Eq. (|4ip correspond to 
the usual Landau-Gilbert description of the Mn spin dy- 
namics induced by the magnetic anisotropy field H. They 



are similar to the results of Refs. [50ll52ll53l and also in- 
clude nonlinear corrections. The first term on the rhs de- 
scribes the Mn spin rotation around the anisotropy field 
H. 7, Eq. (|^51) . gives the renormalization of the gyromag- 
netic ratio by the magnetic exchange interaction. The 
second term describes the Gilbert relaxation of the Mn 
spin towards the direction of H, with Gilbert damping 
coefficient oc a, Ea. (|l^ . The precise values of 7 and 
a in the realistic system depend on the details of the 
valence bandstructure. The relaxation of the hole spin 
component parallel to S{t), s|, depends on the two fields 
h{t) (nonlinear optical excitation) and H(i) (magnetic 
anisotropy) as described in Appendix [D] 

The nonlinear optical excitation gives rise to two ad- 
ditional contributions to the Mn spin equation of motion 
Eq. (|¥T1) . These contributions describe the response of the 
magnetization to the effective magnetic field pulse h(i), 
which is generated by the optical excitation via second 
order processes. The third term on the rhs of Eq. (j4ip 
describes a Mn spin rotation around h{t), whose magni- 
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FIG. 8: (Color online) Light-induced effective magnetic field 
components for pulse durations and Rabi energies as in Fig[T] 
l/rx=21fs, r2 = 10.5fs. Solid line: h^{t), Dashed line: hy{t), 
Dotted line: hz{t). 



tude is proportional to the Gilbert damping coefRcient 
a, Eq. ((44)l . and therefore vanishes if r_L = (hole spin 
conserved) . The last term describes a spin relaxation to- 
wards h(t), determined by the exchange interaction. The 
above light-induced rotation and relaxation vanish after 
the decay of the optical pulse and e-h polarizations. 

Eq. (PT|) predicts an extremely fast response of the spin 
system, determined by the effective magnetic field pulse 
h(i), whose speed is limited only by the pulse duration. 
This result points to a new way of controlling the mag- 
netization during femtosecond time scales, much shorter 
than those accomplished so far. The last two terms on 
the rhs of Eq. (PT|) demonstrate an ultrafast modification 
of the effective magnetic field H. Eq. lHTI) therefore pre- 
dicts two regimes of time evolution, which can be sepa- 
rated due to their very different timescales. The initial 
temporal regime is controlled by the pulsed field h{t) 
and vanishes soon after the optical excitation. The sec- 
ond temporal regime is controlled by the anisotropy field 
H(i), Ea. (|40|) . due to the thermal Fermi sea, which de- 
pends nonlinearly on S(t). 

We now turn to the effective light-induced magnetic 
field pulse h(i). As can be seen from Ea. (|^T|) . in order to 
trigger magnetization dynamics, h(i) must have a com- 
ponent perpendicular to the ground state Mn spin and 
the easy axis. The direction of h{t) is determined by the 
optical transition selection rules and by the interband 
polarizations. Figs [5] and [5] show the three components 
h.x{t), hy(i), and h.z{t) in the case of strong spin dephas- 
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FIG. 9: (Color online) Same as Fig|8] but with rx=0, 
r2=330fs. 



ing (FiglH]) or uncoupled hole spin states (Figl9|). ft is 
clear by comparing the two figures that the direction of 
h{t), which governs the Mn spin trajectory, depends on 
the dephasing. The mosty striking effect of Fj, is the de- 
velopment, on the time scale of the optical excitation, of a 
large component in the direction of optical pulse prop- 
agation. This component is strongly suppressed in the 
case of weak dephasing. As can be seen from Ea. (|4T|) . 
rotates the Mn spin within the x-y plane, in a counter- 
clockwise direction, and leads to spin relaxation towards 
the positive z-axis, i.e. to a magnetization tilt out of the 
x~y plane. On the other hand, for T± — 0, the Gilbert 
damping coefficient vanishes, a = 0, which suppresses 
the light-induced Mn spin rotation. The only effect of 
the photoexcitation is then the relaxation described by 
the last term of Eq. (|4T|) . which causes the Mn spin to 
relax towards the negative y-axis (counter-clockwise ro- 
tation within the x — y plane) since « in this case. 
The above trends are consistent with the full numerical 
calculation of the mean field equations discussed in the 
previous section. In Appendix |E] we derive an equation 
of motion for h(t) that demonstrates its dependence on 
the optical transition selection rules, the hole spin polar- 
ization, the magnetic exchange interaction, and the hole 
spin dephasing. 



VIII. CONCLUSIONS 

In conclusion, we presented a theory of nonlinear spin 
dynamics in ferromagnetic semiconductors and demon- 
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strated the main effects common in all such systems by 
considering a single spin-degenerate electron and hole 
band. Our main result of relevance to experiments is the 
existence of two different temporal regimes of magneti- 
zation dynamics. The initial regime lasts for time scales 
comparable to the optical pulse duration. It is governed 
by magnetization precession and relaxation around an ef- 
fective magnetic field pulse generated by the optical ex- 
citation via second-order coherent nonlinear processes. 
The second temporal regime is governed by the magnetic 
anisotropy due to the thermal carriers. We showed that 
the shape of the magnetization trajectory depends sensi- 
tively on the hole spin dephasing, the magnetic easy axes, 
the transition matrix elements, and the optical pulse in- 
tensity, duration, and helicity. The latter determine the 
magnitude and direction of the ultrafast magnetization 
tilt from the magnetic easy axis in response to the pho- 
tocxcitation of hole spin, which in turn determines the 
importance of nonlinear magnetic effects during many 
picoseconds. To interpret our numerical results, we de- 
rived an effective Landau-Gilbert-like equation of mo- 
tion. This equation shows that an effective magnetic field 
pulse, generated by the interband optical polarizations 
via a Raman-like process, triggers precession and relax- 
ation during femtosecond time scales. This magnetic field 
pulse may be thought of as an ultrafast correction to the 
anisotropy field. The picture of ultrafast magnetization 
dynamics conveyed by our results and, in particular, the 
existence of two temporal regimes is consistent with the 
recent experimental observation of a light-induced mod- 
ification of the magnetic anisotropy field in Ga(Mn)As 
during femtosecond timesales, induced by nonthermal 
high energy carriers A more detailed comparison be- 
tween theory and experiment must be performed in the 
future. Our calculations point out the need for further 
experiments in order to explore the control of the mag- 
netization trajectory by changing the optical excitation 
intensity, duration, central frequency, and helicity. For a 
complete understanding of the magnetization dynamics 
in ferromagnetic semiconductors, the role of fluctuations 
and carrier-spin correlations,^- the bandstructure, and 
disorder must also be considered. 

This work was supported by the EU STREP program 
HYSWITCH. 



APPENDIX A 

In this Appendix we derive the relaxation and dephas- 
ing contributions to the hole spin equations of motion. 
The hole spin is expressed in the basis of hole spin eigen- 
states as 



1 



(^-kmr^^-k^T) - ('l-kmi'l-kmi) , (Al) 



'km+ 



IS 



-km 



i)- (A2) 



From Egs. lT^l) and pp)) we obtain for the corresponding 
hole density matrices after using straightforward algebra 



and the property = 0: 

9t (^-kmJ|^-k™t> 



relax 



■('i-km^'i-kmt) (A3) 



By choosing the z-axis parallel to the Mn spin and noting 
the expression Eq. (jA2p . we obtain Eq. ((23|) that describes 
the dephasing of the hole spin component s^m± perpen- 
dicular to the Mn spin. Similarly we obtain from Eqs. (|19p 
and (Uni) 



relax 



Fj. 



(^-km1>^-kmir) 



"('''-kmt'''-km4^ - kmJI ^ -kmft ) 



which conserve the total hole population 



(A4) 



A^km - (^tkm4^-km4> + {h\,ni,h -^.m^ ■ (A5) 

Using the factorization 

{h\hlhj,h) = {h\h){hlh) ~ {h\h3){hlhi) (A6) 
and the relation (obtained from Eqs. (|Al| ) and (|A5p i 



(A7) 



where cr = ±1, we obtain Eq. ([M|) for the hole spin relax- 
ation. 



APPENDIX B 

In this appendix we derive the polarization dephas- 
ing dtPk,nnaleiax ^siug Eqs.l^ and Af- 
ter straightforward algebra and using the properties 
-^km^kmcr-Z^km = 0, fi^ = 0, we obtain after noting that 
the terms in the summation of Ea. (|19p vanish unless 
(k'm') = (km): 



dt.Pkmnalrelax ~ ^ {h_^^,^^h^k.mi[h^]imaekn) 

F 



^(^-kmo-ekn^-km^^Lkm^)- (Bl) 



We now note the relation between the hole operators with 
spin along the z-axis of pulse propagation or along the 
Mn spin direction: 

hkmit = e'"^^^ cos - /ikmT + e"*"^/^ sin - hkmi , (B2) 

6 
hkmiy = -e'"^/^ sin - h^m^ + e-*"^/^ _ f^^^^ (33) 

where 9{t) and 0(<) are the polar coordinates that define 
the direction of the Mn spin S(t) and the operators /ikmT 
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and humi create hole states with spin along the z-axis. 
Using Eqs. (|B2]) and (|B3]), the property h'^ = 0, and the 
relations / S = e*"^ sinO, Sz/S — cos9 we obtain from 
Eq. ljBip after straightforward algebra 



-r 



5 \n-_kmi"-kmi/l-kmTek«i) 



-r. 



l + Sz/S 



(/l-kmiekni^-kmi^Lkmi) 



■^(^-kmT^kni^-kmi^Lkmt) 



and 



StPkmnilrg;^^ 



(B4) 



I + Sz/S 



+ 2^('i_kmi^-kmT^-kr«iek,ix) 



-r 



1 - Sz/S 



(^-kmiekriJ,/l-kmT'*-kmT) 



2^ ( '^-kmj.eknx/l-kmt^-kmi) 



(B5) 



We now factorize the higher density matrices in the above 
equations as follows: 

{h\h2he) = {h\h2){he) - {h\h){h2e), (B6) 
{h3eh2h\) = {he){h2h\) - {h2e){hh\). (B7) 

We thus obtain that 



I 2 ^ ("-kmi-fkmnT ^ ^km--Pk 
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V k:m+ -Pk 



1 + Sz/S J ^ 

;^ [-PkmnT(l ^ "-krnx) + ^kmnXS^m-] 
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.rn'\ ) 



(B8) 



and 



dtPi 



kmnX I relax 

r± < ^ ("kmX^kmnX ^ S^„j^Pk„i„x) 

(^km-^kmjiX ~ '^'kmX^kmnx) 

^ [-Pkmnx(l — '^kmx) + -PkmriTSkm+] 

[Skm-^icrnnX + ^kmntCl " "k^x)] } ' (^9) 



where n\^m.(j are the spin-polarized hole populations. Af- 
ter expressing the latter by using Ea. (|A7p we obtain af- 
ter some algebra the polarization dephasing contributions 
Eqs. ^ and 



APPENDIX C 

In this appendix we derive an effectve equation of mo- 
tion for the Mn spin by expanding around the adiabatic 
limit. First we use the hole spin decomposition Eq. (PT|) 
into a component parallel S{t) and a component 

perpendicular to S{t) (sjj^). The unadiabatic contribu- 
tion sj^j_ is responsible for triggering the Mn spin preces- 
sion. After substituting Eq. (PT|) into Eq. fTSl) and noting 
that the Mn spin magnitude 5 remains constant in time 
we obtain 



dts^^ = /3cS X s(j^ + Im hk(t) - sl«dtS 



II l-^k|| 



r 



(Cl) 



By projecting out the components parallel and perpen- 
dicular to S and using the relations S ■ dtS = dtS^ /2 = 
and sjj^ • S = we obtain that 



and 



/m /ik||(<) + s^^ -atS (C2) 



Im hk± -r^sjj^ -4||at^C3) 



We now focus on the perpendicular component s']_. We 
note from Eq.® that the motion of S(i) is characterized 
by a precession frequency uj determined by the field H 
and the mean hole spin /3s^: 



5S = w X S , w = /3si - H 



where 



(C4) 



(C5) 



is the mean hole spin. In III-Mn-V semiconductors, this 
precession is much slower than the motion of the hole 
spin s^, which is characterized by the precession energy 
(3cS and the relaxation rate T±. When considering the 
slower Mn spin dynamics, we can then substitute by 
its steady state value in the frame of reference that ro- 
tates with the Mn spin and S(t). The rate of change 
of any vector A as seen by an observer in the rotating 
frame, 9tA|^^j, is related to the corresponding rate of 
change in the inertial frame, 9tA, by 



9* A = dtA[ 



X A. 



We thus obtain in the case of the hole spin 



(C6) 



(C7) 
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Projecting in the direction perpendicular to the Mn spin 
using Eq.([221) we obtain that 

dA^\^ = dtsll^^^^ + H^^^i^, (C8) 

where we used the vector property 

A X (B X C) = (A • C)B - (A-B)C. (C9) 

The rotating frame approximation corresponds to ne- 
glecting dts'll^^^^. 

To obtain an effective equation of motion for the Mn 
spin, we use Eq.® for s'j^ x S to eliminate the hole spin 
s'[, defined by Eqs ^ and from Eq.(IC3|: 

, , s X atS + s x (h X s) 

s'l = S X [s'lxSj = ^— ^ ^.(ClO) 

By summing Eq. (|C3|) over all momenta we thus obtain 
within the rotating frame approximation after using the 
definition Eq.lgl]) that 

(^pcS + /3s| + = {(3cS + S X H 

+/3S X (h X + Tj_S X (SxH- ^tS) . (Cll) 

The above nonlinear equation of motion is governed by 
the magnetic field H, the effective magnetic field h deter- 
mined by the interband e-h polarizations, and the hole 
spin relaxation (Fji) and dephasing (r_L) rates. Eq. (|Cll|) 
can be transformed into an effective time-dependent 
Landau-Gilbert-like equation of motion by taking the 
cross product of both sides with S and using the vector 
property Eq. (jC9p and the property S ■ dtS — 0: 

(^(3cS + /5s|f + dtSxS^ {(3cS + S x (h x s) 
+l3hxS + T^S xU-T^dtS. (C12) 

Substituting the above expression into Eq. (|Clip we ob- 
tain Eq. fiT]) . 

APPENDIX D 

In this appendix we derive the relaxation of the hole 
spin component s|| parallel to the Mn spin S(i), deter- 
mined by Eq. (|C2[) . We eliminate s'[ from Eq. (IC2p by 
using the property • dtS — H • dtS//3, obtained after 
using the properties {SxdtS)-dtS = 0, S-dtS = dtS"^ = 0, 
and some algebra. We also eliminate 9tS from Eq. (jC2[) 
by projecting both sides of Eq. dlTI) with H and using the 
vector property A • (B x C) = B • (C x A) and the re- 
lations H2 = Hj -I- 2 and h • H = h\\H\\ • H_l. 
By summing both sides of Eq. (jC2[) over all momenta we 



finally obtain the equation of motion 

9.4 + r,, (s\ + m.) + r^Hi+/3H^ .h^ 

(/3c5 + /34 + i/||j +Tl 

ri(Hxh)|| -(/3c5' + £f||)Hi-h^ 
= h\i + -2 . (Dl) 

(pcS + f3sl + HA +Tl 



APPENDIX E 

In this appendix we derive the equation of motion 
that determines the time evolution of the effective ma- 
gentic field pulse h{t) = £{t)h{t) that governs the in- 
tial femtosecond magnetization re-orientation and relax- 
ation. Using the equations of motion of the interband 
polarizations, Eqs. (PT|) . ([5^. ([M]). and ([33]), we obtain 
the equation of motion 

z9thk = r!khk + ^ Ak X hk + + £it)Sk, (El) 

where fik and Ak were defined by Eqs. ([35]) and ([35]) . 

PkW-A^;^kTW+'^-^kiW (E2) 
satisfies the equation of motion 

idtp^ - rikPk + - 5(i)A/-k, (E3) 

where 

AAk=(|/.+ p + |M-n [l-(iV^ + /2] 

-(Im+P-Im-P)G^L-4.), (E4) 

and 

5k. - (Im+P + Im-H 4. + * (Im+P - Im-P) si^y 

S^y = (Im+P + Im-H s'iy ^ (Im+P - Im-H 4. 
+2Re{fi+fi*_)sly - 2/m(^+/il)4^, 

5k.--(|/i+p-|M-P) [l-{Nt + K)/2] 

+ (Im+P + Im-P)(4.-sL)- (E5) 

We note from the above equations that, in the absence of 
spin polarization in the ground state, 8 = 8^=8*^ = 0, 
h(i) points along the z-direction of optical field prop- 
agation. On the other hand, in the presence of spin- 
polarized holes in the ground state as in ferromagnetic 
semiconductors, h develops additional components de- 
termined by 8'', the Mn spin and exchange interaction 
Ak, and by the magnetic anisotropy. 
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